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Abstract 



We study the two-dimensional fully frustrated XY (FFXY) model and two 
related models, a discretization of the Landau-Ginzburg- Wilson Hamiltonian 
for the critical modes of the FFXY model and a coupled Ising-XY model, 
by means of Monte Carlo simulations on square lattices L 2 , L < 10 3 . We 
show that their phase diagram is characterized by two very close chiral and 
spin transitions, at T c h > T sp respectively, of the Ising and Kosterlitz-Thouless 
type. At T c h the Ising regime sets in only after a preasymptotic regime, which 
appears universal to some extent. The approach is nonmonotonic for most 
observables, with a wide region controlled by an effective exponent ~ 0.8. 

PACS Numbers: 74.81.Fa, 64.60.-i, 64.60.Fr, 75.10.Hk 
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The nature of the phase transitions in frustrated systems is of great interest in statistical 
physics. The two-dimensional fully frustrated XY (FFXY) model [1] is defined by the 
Hamiltonian 

Hffxy = -J^2cos(9 x -9 y + A xy ), (1) 

where the sum is over all nearest-neighbor pairs of a square or triangular lattice, and A xy 
satisfy the constraint A xy = tt around any plaquette. It is experimentally relevant for 
Josephson-j unction arrays in a magnetic field [2]. The angle variable 9 X corresponds to the 
phase of the superconducting order parameter on each superconducting grain, and A xy is 
the vector potential of a perpendicular magnetic field corresponding to half a flux quantum 
per plaquette. 

The ground state of the FFXY model presents an enlarged 0(2)£g>Z 2 degeneracy. The 
additional Z 2 degeneracy is related to the breaking of the chiral symmetry [1]. For each 
plaquette II we define the chirality 

C n = MOx-Oy + Axy), (2) 

(xy)eU 

where n is the dual-lattice site at the center of II. On the square lattice the staggered 
magnetization Mc = ^2 n (— l) ni+U2 C n /V defines an order parameter, which competes with 
the spin modes to determine the phase diagram of the FFXY model. The critical behavior 
has been much investigated during the last few decades (see, e.g., [3-21]) using Monte 
Carlo (MC) simulations, real-space renormalization-group (RG) techniques, field-theoretical 
methods, etc. Moreover, several related models have also been considered: for instance the 
fractional-charge Coulomb gas, coupled Ising-XY models, coupled XY models. 

In spite of all this work, the phase diagram and critical behavior of the FFXY model 
are still rather controversial. Most recent MC simulations favor the existence of two very 
close transitions [7,8,10,12,13,15,16,19]. The most likely interpretation is that the higher- 
temperature transition is characterized by the onset of chiral long-range order, while spins 
remain disordered. The lower-temperature transition is associated with the breaking of the 
continuous symmetry and is followed by a low-temperature phase in which spin quasi-long- 
range order coexists with chiral long-range order. The chiral transition is expected to be 
in the Ising universality class, due to the scalar nature of the chiral order parameter. The 
second one should be a Kosterlitz-Thouless (KT) transition. This scenario is also supported 
by arguments based on a kink-antikink unbinding picture [17,21]. The Ising nature of the 
chiral transition has not been satisfactorily supported by numerical simulations so far. Most 
MC simulations [4-9,11,12,14,15,19,20] have found that the behavior of the chiral modes — 
both in the FFXY and in related models — is not consistent with an Ising transition. For 
example, most finite-size scaling (FSS) analyses have obtained v 0.8, instead of the Ising 
value v — 1. There are several possible explanations. One possibility is that, even if the 
chiral order parameter is a scalar, the chiral transition belongs to a universality class that 
is not the Ising one. After all, the estimate v ~ 0.8 appears to be somewhat universal, the 
same value being obtained in several different models. A second one is that the Ising regime 
sets in only on large lattices [10]: the observed behavior is only an intermediate crossover. A 
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FIG. 1. Sketches of the phase diagrams of the 4> 4 model for U = 1 (above) and of the IsXY 
model (below). The behavior in the circled region where the transition lines meet is unclear. 



third possibility is that the apparent two transitions are a finite-size effect. On the contrary, 
spin and chirality order at the same temperature. In this case the critical spin and chiral 
modes couple at criticality and give rise to a qualitatively new critical behavior, in which 
chiral modes do not behave as spins in the Ising model. 

In this paper we consider the square-lattice FFXY model (implemented by alternating 
vertical lines with ferromagnetic and antiferromagnetic couplings) and two related lattice 
models. The first one is a 4 model defined by the Hamiltonian 

= - j ■ <kv + E + u^l - !) 2 ] 

(xy),i i,x 

+2(U + D)J2<Plx<Pl,x, (3) 

X 

where <f>i tX , (p2,x are real two-component variables denned on the sites x of a square lattice, 
0? = (f)i ■ (pi, and J,U,D > 0. Hamiltonian Ti^ is expected to describe the critical modes of 
the FFXY model. It corresponds to a straightforward lattice discretization of the Landau- 
Ginzburg- Wilson theory obtained by applying a Hubbard- Stratonovich transformation to 
the FFXY model, see, e.g., Ref. [5]. Here we set U — 1; all results reported below refer to 
this case. The symmetry <pi <-> 2 is the analogue of the Z 2 chiral symmetry of the FFXY 
model. The chiral order parameter is 

C x = <Pl, x - <Pl, x - (4) 

When D = model (3) becomes 0(4) symmetric; it does not have any transition at finite 
temperature, but its correlation length diverges exponentially for J — > oo, see, e.g., Ref. [22]. 
We also consider the coupled Ising-XY (IsXY) model [4] 

T^isXY = - ^ — (1 + cr x a y ) s x ■ s y + Ca x a y , (5) 

(xy) 



TABLE I. Comparison of estimates of r) obtained from % s , and from R s and T using the 
relations valid in the XY model, such as Eq. (6), in the low-temperature phase. 

model J V V from R s r\ from T 

FFXY 2A 0.1480(5) 0.1479(4) 0.14779(11) 

2.3 0.1750(5) 0.1752(4) 0.1758(3) 

2.26 0.2023(11) 0.2015(7) 0.2021(5) 

(j>, D = 1/2 1.50 0.1341(6) 0.1346(4) 0.13425(12) 

1.48 0.1675(16) 0.1672(10) 0.1670(4) 

<t>,D = 99 1.6005 0.120(3) 0.1214(6) 0.1202(5) 

IsXY, C = 1.52 0.1817(7) 0.1821(4) 0.18190(13) 



where a x = ±1, and the two-component spins s x satisfy s x ■ s x = 1. Here s x and a x 
correspond to spin and chiral variables, respectively. Note that, by taking the limit U — > oo 
and then D — > oo in model (3), one recovers the IsXY model for C = 0. 

We performed MC simulations on L x L lattices with periodic boundary conditions, 
and sizes up to L = O(10 3 ). We used mixtures of Metropolis and overrelaxation updating 
algorithms, as proposed in Ref. [23]. In this paper we present the main results; details will 
be reported elsewhere. We find two very close Ising and KT transitions in the FFXY model, 
and in the 4 and IsXY models for extended ranges of the parameters D > and C. No 
evidences of unique continuous transitions are found. For sufficiently large D and — C, i.e. 
D > 50 and C < —5, we find instead a single first-order transition. The phase diagrams of 
the 4 and IsXY models are shown in Fig. 1 [24]. Little is known about the region where 
the two continuous transition lines turn into a single first-order one. Another interesting 
feature emerges from our numerical results: the asymptotic critical behavior at the chiral 
transition is observed only after a peculiar nonmonotonic crossover regime, which appears 
universal to some extent. 

In the FFXY model the square lattice can be divided into four sublattices, so that the 
four sites of every plaquette belong to different sublattices. The ground state is translation 
invariant within these sublattices. We define the spin correlation function G s (x) as the 
correlation between two spins in the same L/2 x L/2 sublattice. Then we consider the 
staggered chirality correlation function G c (x). In the <p A model the spin and chiral correlation 
functions are defined by G s (x) = 0i,o'0«,:r) and G c (x) = (Co C x ) c . In the IsXY model we 
have G s (x) = (s • s x ) and G c (x) = (o" o" x ) c . From G s and G c we define the susceptibilities 
Xs and Xc and second-moment correlation lengths £ s and £ c . We also consider R s = £ S /L, 
R c = £ C /L, the spin and chiral Binder parameters B s and B c , and the helicity modulus T. 

We first show that the low-temperature phase of these models is characterized by the 
breaking of the Z 2 chiral symmetry and by a spin quasi-long-range order analogous to the 
one of the standard XY model. Indeed, chiral modes are magnetized and, in the large-L 
limit, the exponent r\ (computed using Xs ~ L 2 ~ v ), R s = £ S /L, and T satisfy the universal 
relations that hold among the corresponding quantities in the low-temperature phase of the 
XY model on a L x L square lattice with periodic boundary conditions. For example [25,26] 

27TT] 7] 2 E„=-oo expt-Tm 2 /??) 

for < r] < 1/4. As shown by the results of Table I, the agreement is very good and provides 
a conclusive evidence that the low-temperature phase of the FFXY and related models is 
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FIG. 2. R s = i s /L at fixed R c = £ C /L vs 1/L. 



controlled by the same line of Gaussian fixed points that is relevant for the XY model. 

To determine the number of transitions, we perform a FSS analysis using the method 
proposed in Ref. [27], see also [28]. Instead of computing the various quantities at fixed 
Hamiltonian parameters, we compute them at a fixed value of a chiral — this guarantees that 
we are observing the chiral transition — RG invariant quantity; in our specific case we choose 
R c = £ C /L. This method has the advantage of not requiring a precise estimate of the critical 
temperature. We fix R c = R~i s where Rj s = 0.905048... is the corresponding Ising value [29]. 
Note that such a choice does not represent a bias in favor of an Ising transition. R c can be 
fixed to any value R c ^: the conclusions are independent of R c ,6x- If the transition is unique, 
also the spin correlation length £ s diverges, and thus R s = (, S /L at fixed R c is expected to 
converge to a nonzero value. If there are two transitions and T c h > T sp , £ s remains finite, so 
that R s ~ 1/L for L — > oo. Results for R s are shown in Fig. 2. They appear to decrease with 
increasing L, without showing any hint at a convergence to a nonzero value. This shows 
that the spin correlation length ^ at the chiral transition is finite, though quite large. For 
example, for L — ► oo we have £i c) = 118(1) for the FFXY model, > 400 for the 4 model 
with D = 1/2, = 52.7(4) for the IsXY model at C = 0. Additional evidence is provided 
by the helicity modulus. At the chiral transition, it appears to vanish in the large- L limit, 
consistently with the fact that spin variables are disordered. 

As a consequence of the large value of £s , the eventual Ising regime at the chiral tran- 
sition may set in only when L ^> £g [10]. This fact is substantially confirmed by our FSS 
analysis. For instance, if, for any invariant quantity S, we define an effective exponent u e g(L) 
as 



l/u eS (L) = (\ndS/dJ\ 2L -\ndS/dJ\ L ) /\n2, (7) 

we find that u^(L) shows first a plateau at u e s(L) ~ 0.8 — this is consistent with previous 
works on small lattices — and only for L > ^ appears to approach the Ising value. Similar 
crossover effects are observed in other quantities. For example, the chiral Binder parameter 
B c is systematically higher than the Ising value [29] Bi s = 1.167923(5). However, as L 
increases, B c starts to decrease and, if L is large enough, its behavior is compatible with the 
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FIG. 3. AB C = B C - Bi a at fixed R c = £ C /L vs L" 7 / 4 . 



expected one [30], B c = B\ s + aL VlB ~ 2 = Bi s + aL~ 7 ' 4: , see, e.g., Fig. 3. This confirms the 
Ising nature of the chiral transition. 

The spin transition, whenever continuous, is expected to be a KT one. We verify that 
at the KT transition R s and T behave as expected, for example [25] 



T = 0.63650817819... + 



0.318899454... 
InL + c 



+ ... 



(8) 



In most cases, including the FFXY model, the chiral and the spin transitions are very 
close. If 8 = (J sp — Jch)/Jch, J S p and J c h being the location of the two transitions, we find 
5 = 0.0159(2) for the FFXY model [J ch = 2.2063(1) and J sp = 2.2415(5)], 5 = 0.0025(2) in 
the 4 model for D — 1/2, 5 — 0.0167(7) in the IsXY model with C = 0. 

Our FSS analysis definitely shows that the chiral transition, when it is continuous, be- 
longs to the Ising universality class. However, the Ising critical regime is reached only after 
a crossover region in which effective exponents and RG invariant quantities show a behav- 
ior that is surprisingly similar in the FFXY model, in the 4 model with < D < 2, 
and in the IsXY model for —0.5 < C < 0.2. Fig. 4 shows the FSS curves of R s = £ S /L 
at fixed R c = £ c / L plotted versus a rescaled lattice size L r = L/l, where I is a rescal- 
ing factor that is chosen judiciously for each different model. The optimal data collapse 
is obtained for Z/Zffxy — 4.6,3.2,1.8,1.1 for the 4 model at D — 1/3,1/2,1,2 and 
l/l FFXY = 2.6, 1.1, 0.45, 0.09 for the IsXY model at C = -0.5, -0.2, 0, 0.2. Of course, for L 
large £ s — > ^ and therefore, given two different models, h/h should correspond to the ratio 
of the corresponding spin correlation lengths at the chiral transition. Thus, the scaling we 
observe implies that R s is an approximately universal function of L/t^ . Fig. 5 shows the 
effective exponent l/^ e ff defined by Eq. (7), as obtained from R c , B c , and R s , respectively, 
by using the same rescaling factors I as determined from R s . In all cases we observe an 
approximate collapse of the data. Note that there is a rather extended region, L r < 1, i.e. 
L < !;i c \ in which z/ e g computed from the chiral variables R c and B c (resp. R s ) is approx- 
imately 0.8 (resp. 0.9). This preasymptotic behavior explains previous estimates v w 0.8 
of the chiral exponent, see, e.g., [4-9,11,12,14,15,19,20]. Actually, the data of l/v e & from 
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FIG. 4. i? s at fixed R c versus the rescaled lattice size L r = L/l. We set I = £^ ~ 118 for the 
FFXY model. 



R c and B c appear to undershoot the Ising value, which is then approached from below. A 
data collapse is also observed for B c itself, the helicity modulus, and the effective exponent 
r) e ff obtained from Xc- In conclusion, the crossover to the asymptotic Ising behavior presents 
interesting features and appears universal to some extent. 

The origin of the apparent scaling behavior is not clear. It might reflect the nearby 
presence of a multicritical point, which could be observed only by performing a further fine 
tuning of the Hamiltonian parameters, where chiral and spin modes become critical at the 
same time. The identification of this multicritical point remains an open issue. 
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